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ABSTRACT 

We report the results of the first self-consistent three-dimensional adaptive mesh refinement magne- 
tohydrodynamical simulations of Population III star formation including the Biermann Battery effect. 
We find that the Population III stars formed including this effect are both qualitatively and quanti- 
tatively similar to those from hydrodynamics-only (non-MHD) cosmological simulations. We observe 
peak magnetic fields of ~ I0~^ G in the center of our star-forming halo at z ~ 17.55. The magnetic 
fields created by the Biermann Battery effect are predominantly formed early in the evolution of the 
primordial halo at low density and large spatial scales, and then grow through compression and by 
shear flows. The fields seen in this calculation are never large enough to be dynamically important 
(with /3 > 10^^ at all times), and should be considered the minimum possible fields in existence during 
Population III star formation, and may be seed fields for the stellar dynamo or the magnetorotational 
instability at higher densities and smaller spatial scales. 

Subject headings: cosmology: theory - galaxies: high-redshift ~ magnetohydrodynamics - stars: for- 
mation 



1. INTRODUCTION 

The nature of the first generation of stars, as well 
as their influence on later structure formation, is a 
fundamental problem in cosmology. A great deal 
of theoretic al progress has b een m a de (see recen t 
revi ews by iBromm fc LarsonI (|2004[ ). iGloveil (|2005D . 
and lCiardi fc Ferraral (|2005f ). In the past decade, cosmo- 
logical hydrodynamic simulations of Population III star 
formation have achieved great success, and significantly 
different numerical methods have produced results that 
agree quite well (JAbcl ct a l. 2002; Bromm et al. 2002; 
lYoshida et al.ll200ilO'Shea fc Norman 2007). These cal- 
culations have given a reasonably clear picture of the for- 
mation process of Population III stars, and have provided 
some constraints on many of their important properties. 

These calculations, while useful, largely ignore an im- 
portant issue; the relevance of magnetic fields in Pop- 
ulation III star formation. Magnetic fields are widely 
observed in our galaxy, in other galaxies, and in galaxy 
clusters, and the origin of these fields is one of the most 
fundamental and challenging problems in astrophysics 
(jCarilh & Tavloill2002t IWidrowll2002D . One possibihty 
is that magnetic fields are created and amplified in the 
first generation of stars and are spread throughout the 
IGM when these stars explode, providing seed fields for 
later generations of stars and for further amplification 
by dynamo effects. If some seed magnetic fields exist be- 
fore Population III stars form, they may help to remove 
angular momentum from the star- forming clouds, signif- 
icantly cha nging the ultimate mass range of these stars 
(jPudritz fc Silk 1989; Davics & Widrow 2000). 

Several groups have examined the importance of mag- 
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netic fields on the evolution of Population III pro- 
tostellar disks us ing analytic o r semi-analytic mod- 
els. T hese include ISilk fc Laneeil ((200 6) . Tan & McKei 
()2004D and lTan fc BlackmanI (|2004i rwho model (among 
other aspects of primordial star formation and evolu- 
tion) dynamos i n primordial accretion disks. Othe r au- 
thors, including 'Flo wer fc Pineau des ForetsI ()2003[ ) and 
[Maki & Susa (2007), use one-zone calculations to exam- 
ine the collapse of the primordial star-forming cloud and 
the assumptions of flux-freezing. While useful, these 
models do not self-consistently include the effects of both 
magnetic fields and cosmological structure formation. 

In this Letter, we report the results of the first mag- 
netohydrodynamic simulations of Population I II star for- 
matio n including the Biermann Battery effect (jBiermannI 
Il950( l within the context of cosmological structure for- 
mation. We describe our numerical methods and the 
simulation used in Section [21 present our key results in 
Section [3l and discuss and summarize our results in Sec- 
tion H 

2. METHODOLOGY 

The simulation discussed in this Letter was performed 
using the a daptive mesh N-bod y plus hydrody namics 
code Enzo (jO'Shea et all 12004 iNorman et aLl 120071 ). 
which has recently been extended to including the equa- 
tions of ideal magnetohydrodynamics. The MHD solver 
is a h igh-order God un ov-type finite-volum e numerical 
solver (iLi fc Lill2003l : 120051 : III et al.ll2008fl . This solver 
was re cently successful ly used to study magnetic jets 
(jLi et a l. 2006; Nakamuxairanr2006, 2007). To main- 
tain the divergence- free magnetic field condition, we used 
a con strained transport (CT) scheme (jBalsara fc Spiceil 
Il999( l as well as a modified div ergence-f r ee rec onstruc- 
tion scheme original proposed bv'Ba lsaral (|200i r) that in- 
cludes second-order-accurate divergence-free restriction 
and prolon gation for magnetic fields in the adaptive mesh 
hierarchy () Collins fc NormanI l2008l ) . These equations 
are solved in cosmological calculations in the comov- 



ing frame, as described in IXu et al] (|2008f ). To track 
the pressure more accurately in hypersonic regions, we 
have im plemented the modified entropy equation given in 
iRvu et al. ( 1993) and the internal energy equation given 
in lBrvan et al.1 (| 19951 ) in our code. 

To study the effect of the Biermann battery during 
the formation of Population III sta rs, the battery term 
is added to the induction equation (|Kulsrudll2005[ ): 

— - = V X (v X B + 2 
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e 1 + X P^ 
where c is speed of light, pe is the electron pressure, rie 
is the electron number density, e is the electron charge, 
m is the average mass per particle and x is the ionization 
fraction. This term is non-zero in regions where curved 
shocks are formed. The second equation assumes that 
m and x are unchanged in space, which i s a reasonable 
appro ximation locally. As pointed out by iKulsrud et al .l 
(|1997t ). the ionization fraction enters into the induction 
equation through 1 + Xj so even a very small ionization 
fraction is enough to generate magnetic fields. In our 
simulation, we took m = 1.2 mn, where mn is the mass 
of a neutral hydrogen atom. Note that the electron frac- 
tion is <C 1 at all times in this star formation scenario, 
so we set x = rather than calculating the electron frac- 
tion in every cell. Given the approximation above, this 
introduces negligible error. The battery term is added to 
the simulation through the EMF term in the constrained 
transport algorithm, ensuring that no divergence of mag- 
netic fields is introduced to the system. 

Since the electron fraction declines rapidly at the halo 
core, ambipolar diffusion may be important in star for- 
mation simulation (Abel et al. 2002). To address this is- 
sue, we compared the ambipolar diffusion timescale with 
the dynamic timescale in o ur simulation, as was done 
by (jO'Shea fc Normaiill2007f ). The ratio is always larger 
than 10 over the density range of interest, suggesting that 
ambipolar diffusion does not play a significant role up to 
the densities that we study. 

Our simulation setup i s similar to that described in 
10 'Shea fc NormanI (|2007f l. We use a N-body plus hy- 
drodynamics simulation with a comoving box size of 0.3 
h~^ Mpc, initiali zed at z = 99 with an Eisenstein fc Hu 
power spectrum ([Eisenstein fc IIulll999t) . and with cos- 
mological parameters fit = 0.04, ilm = 0.3, flA = 0.7, 
h — 0.7, ag, — 0.9, and ris — 1.0. The simulation was 
first run with a 128"^ root grid and three levels of AMR 
to z — 15, where the most massive halo in the volume 
was found. The simulation was then re-centered on this 
halo, using a 128'^ root grid and three levels of static 
nested grids, giving a dark matter and baryon resolution 
of 1.81 h"^ Mq and 0.28 h"^ Mq, respectively, and an 
initial comoving spatial resolution of 293 h~^ pc. The 
simulation was started at z = 99, initialized with zero 
magnetic fields throughout, and allowed to run with a 
maximum of 22 levels of adaptive mesh refinement un- 
til the collapse of the most massive halo, at z c± 17.55, 
using the n ine-species mol e cular chemistry of lAbel et al.l 
([1997) and lAnninos et"a L (199'<^, bu t mo dified for high 
densities as described in .Abel et al.l ([2002) . Further de- 
tails of the simulation setup and physics can be found in 



10 'Shea fc NormanI ([200l . 

3. RESULTS 

Over the course of the simulation, magnetic fields are 
produced at various physical scales by the Biermann Bat- 
tery process. Figured) shows the evolution of spherically- 
averaged radial profiles of the baryon number density and 
magnetic field strength of the primordial star-forming 
halo at several times during the halo evolution. Since the 
magnetic fields are very small (/3 = Pthermai / Pmagnetic 
is greater than 10^^ at all times), the gas collapse at the 
halo core is almost the same as the hydrodyn amic {B = 
0) case described bv lO'Shea fc NormanI ([20071 ). The pro- 
files start at z = 40, where the magnetic fields generated 
by the Biermann Battery are of order 1 ~^^ G, which 
is consistent with theoretical p redictions ([Pudritz fc Silkl 
119891 : iDavies fc Widrowl 12000 ') . In plasmas with a large 
beta {'S> 1) and large mass-to-flux ratio, magnetic fields 
follow the gas falling into the core passively, so after be- 
ing generated at low densities, the magnetic fields are 
primarily amplified by being carried along with collaps- 
ing baryons. This is shown more clearly in Figure [2l and 
will be discussed in more detail later. At low density (up 
to nf, ~ 2 cm^-^, corresponding to the evolution of the gas 
up to z ~ 20), the Biermann Battery is effective at cre- 
ating magnetic fields. The central magnetic field scales 
with density as \B\ oc p, suggesting that the magnetic 
field is enhanced by both the Biermann Battery effect 
and by gas collapse. At higher densities (after z ~ 20), 
the magnetic field at the center of the halo is amplified 
as \B\ oc p^^"^, as expected in spherical collapse. At later 
times, when the gas in the center of the halo collapses, the 
magnetic field in the center of the halo rises to ^ 10^^ G 
as the central number density grows to Ub ~ 10^" cm~^. 

It is useful to examine the history of magnetic field am- 
plification inside the halo in more detail. This is shown 
in Figure [21 which shows the mass- and volume- weighted 
magnetic fields inside the halo virial radius as a function 
of redshift (top panel), and also shows the mass- weighted 
mean magnetic field strength in the inner 100 Mq of 
gas as a function of the mass-weighted baryon number 
density of gas in the same Lagrangian region (bottom 
panel). From these two plots, the two stages of mag- 
netic field amplification can be observed. From z = 40 to 
z = 20, the magnetic fields slowly increase to 10~^^ G. At 
this stage, the mass- weighted and volume- weighted fields 
are quite similar, because the magnetic fields are created 
and amplified by the Biermann Battery at large spatial 
scales, and the fields ar e distributed uni f ormly by merg- 
ers and shear flows (see iRoettiger et all ([19990 ). During 
the second stage of evolution, after z = 20, the volume- 
weighted and mass-weighted fields evolve in very differ- 
ent ways, with the mass- weighted magnetic field growing 
rapidly along with the increasing baryon density in the 
halo cores. 

The bottom panel of Figure [21 shows that at late times 
(high densities) , the magnetic field in the core of the halo 
increases as \B\ ex p'^^^, or solely from field amplification 
due to magnetic fields frozen into spherically-contracting 
gas. The collapse timescale at the center of the halo is 
proportional to the cooling time, which is quite short 
at high densities, making any significant magnetic field 
creation by the Biermann Battery mechanism difficult 
in the halo center. The baryon density in the halo core 
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Fig. 1. — Evolution of spherically-averaged, mass-weighted 
radial profiles of baryon number density (top) and magnetic 
field strength (bottom) of the Population III star-forming halo. 
Lines correspond to (from bottom to top in each panel) z = 
40, 30, 25, 20, 19, 18, 17.61, 17.55. 

grows about 9 orders of magnitude (from 1 to 10^" cm~'^ 
between z — 20 and z = 17.55, when the simulation is 
terminated) , and the magnetic fields in the core increase 
by about 5 orders of magnitude (from 10~^^ to 10~^ 
G). Though the core region contains less than one per- 
cent of the halo gas, the mass- weighted magnetic field in 
the halo increases by approximately four orders of mag- 
nitude, from ~ 10^^^ to 10^^^ G, between z=20 and 
z=17.55. 

To show where magnetic field generation via the Bier- 
mann Battery takes place, and to provide a more quan- 
titative understanding of its effect in different stages of 
halo evolution, we show in Figure [3] the two-dimensional 
distribution of the instantaneous magnetic field genera- 
tion rate via the Biermann Battery versus baryon over- 
density at four simulation outputs from z = 40 to 
z — 17.55. At early times, very little gas is at high over- 
densities, and the magnetic field generation rate is below 
10~^^ G/s. The total time taken during the simulation, 
from z — 99 to z = 17.55, is approximately 6.2 x 10^^ 
seconds. If the magnetic fields are only created by the 
Biermann Battery without any amplification from shear 
flows or collapse, the magnetic fields would be at most 
10-19 G at z ~ 18. By z = 40 (~ 1.5 x IQi^ seconds, 
or 5 X 10^ years, after the beginning of the calculation), 
fields reach a strength of 10^ ^^ G, suggesting that some 
amplification is taking place through gas collapse. While 
the rate of magnetic field generation via the Biermann 
Battery can actually be much higher at high overden- 
sity (as seen in the bottom right panel, at z = 17.55), 
the rate never exceeds 10^^° G/s (~ 3 x 10"^^ G/yr). 
Given that the collapse of the core, from a peak baryon 
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Fig. 2. — Top panel: mass-weighted (dashed line) and volume- 
weighted (solid line) mean magnetic field strength inside the virial 
radius of the collapsing halo as a function of redshift. Bottom 
panel: mass-weighted mean magnetic field strength as a function 
of baryon number density inside the core of the primordial star 
forming halo at the final simulation output (solid line; Menc — 
100 Mq). The dashed line shows |-B| oc p^''^ (as expected by flux 
freezing and spherical collapse). 
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Fig. 3. — Mass- weighted two-dimensional distributions of the 
instantaneous rate of magnetic field generation via the Biermann 
battery vs. baryon overdensity, at 2; = 40 (top left), z = 30 (top 
right), z = 20 (bottom left) and z = 17.55 (bottom right). Each 
panel includes all gas within the virial radius at a given redshift. 
Contour lines are logarithmic, and show relative values between 
the highest and lowest non-zero values in each panel. 

density of Uf, '^ 1 cm"^ to 10^" cm"^, takes approxi- 
mately 10"'^'* seconds (~ 3 x 10^ years), magnetic fields 
larger than 10^ ^^ G cannot be made with the Biermann 
Battery during this time, strongly suggesting that the 
battery effect makes little contribution to magnetic field 
amplification during the final collapse. 

4. DISCUSSION AND CONCLUSIONS 

In this Letter, we have shown results from the first cos- 
mological magnetohydrodynamical simulation of Popula- 
tion III star formation including the Biermann Battery 
effect. This effect is one of the most robus t methods 
of generating magnetic fields in the Universe (jBiermannl 



Il950f ). and thus provides useful constraints on the mini- 
mum magnetic field expected in situations where Popu- 
lation Iff star formation will take place. 

We find from our simulation that small magnetic fields 
are primarily generated via the Biermann Battery at rel- 
atively low overdensities, and are then amplified to values 
of nearly 10"^ G at the center of the cosmological halo 
via gravitational collapse. While significant, this mag- 
netic field is still quite small - the plasma f3 is never 
smaller than 10^^ at any point during the simulation. 
This suggests strongly that the magnetic fields do not 
play a significant dynamical role up to n^ ^^ fO^*' cm~^, 
when the simulation is terminated. As a result, the evo- 
lution of the primordial gas within the cosmological halo 
is very close to the hydr odynamic results obtained by 
lO'Shea fc NormanI (|2007D . At later times, however, the 
magnetic fields may be amplified to dynamically relevant 
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values via a process such as the magnetorotational insta- 
bility or a dvnamo (|Sillt fc Li 
I2OOI Wa^k Blackmanll2004D . 

The simulation of further evolution of our evolving pro- 
tostellar core is currently beyond the capability of our nu- 
merical tools. However, there are multiple avenues that 
we can follow to more thoroughly explore the relevance of 
magnetic fields in Population III star formation. We can 
begin our calculation with seed fields having s trengths 
based on limits from observation and theory (jWidrowl 
l2002f) . We can examine the effects of magnetic fields 
at higher densities by implementing more chemical an d 
cooling processes (as discussed by I Turk et all I2008D . 
Both of these projects are underway, and we will report 
on the results in an upcoming paper. 



To summarize, we have performed the first calculations 
that incorporate the Biermann Battery in cosmological 
magnetohydrodynamic simulations of Population III star 
formation. Our key results are as follows: 

1. From an initial state with no magnetic fields, a 
combination of the Biermann Battery and compressional 
amplification can result in fields with strengths of \B\ ~ 
10"^ G ntuB — W^° cm~'^ at the center of a cosmological 
halo where a Population III star will form. 

2. The Biermann Battery creates fields predominantly 
at low density (n^ < 10 cm~'^) and large spatial scales 
in Population III star-forming halos. 

3. The magnetic fields created by the Biermann Bat- 
tery are dynamically unimportant at all densities below 
ns ~ 10^° cm-3 -13 = Pth/Ps > 10^^ at aU times during 
the evolution of the halo. 
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